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ABSTRACT 

Context. Clouds containing molecular dark matter in quantities relevant for star formation may exist in minihaloes of the type of cold dark 
matter included in many cosmological simulations or in the regions of some galaxies extending far beyond their currently known boundaries. 
Aims. We have systematically explored parameter space to identify conditions under which plane-parallel clouds contain sufficient column 
densities of molecular dark matter that they could be reservoirs for future star formation. Such clouds would be undetected or at least appear 
by current observational criteria to be uninteresting from the perspective of star formation. 

Methods. We use a time-dependent photon-dominated region code to produce theoretical models of the chemistry and emission arising in 
clouds for our chosen region of parameter space. We then select a subset of model clouds with levels of emission that are low enough to be 
undetectable or at least overlooked by current surveys. 

Results. The existence of significant column densities of cold molecular dark matter requires that the background radiation field be several or 
more orders of magnitude weaker than that in the solar neighbourhood. Lower turbulent velocities and cosmic ray induced ionization rates than 
typically associated with molecular material within a kiloparsec of the Sun are also required for the molecular matter to be dark. 
Conclusions. We find that there is a large region within the parameter space that results in clouds that might contain a significant mass of 
molecular gas whilst remaining effectively undetectable or at least not particularly noticeable in surveys. We note briefly conditions under 
which molecular dark matter may contain a dynamically interesting mass. 

Key words. ISM: clouds - ISM: molecules - galaxies: ISM - galaxies: kinematics and dynamics - stars: formation 



1. Introduction 

The mapping of external galaxies in the 21cm Hi (e.g. 
iNoordermeer et alJl2005l IWalter et alJl2005l) and the 2.6 mm 
CO fe.g. iHelfer et alJl2003t Lundgren et aljEooH) lines pro- 
vides information about the distributions of atomic and molec- 
ular gas in those galaxies. Many studies of these distributions 
are undertaken in the context of star formation. For example, 
one commonly addressed question concerns whether a column 
density of gas in excess of a local critical value leads to star 
formation ( Kennicutt 1989). 

Observations of far infrared fine structure lines indicate that 
in metal-poor regions, like those found in some dwarf galax- 
ies, 2.6 mm CO observations alone lead to a significa nt un- 
derestimate of the mass in molecular gas (Mad den et alJll997t 
Madden 2002). This is consistent with the results of multi- 
line CO studies showing that the 2.6 mm integrated antenna 
temperature to molecular hydrogen conversion factor in such 
regions is much larger than is appropriate for Galactic giant 
molecular clouds, and also with the known difficulty in detect- 
ing CO emissions in some dwarf and compa ct galaxies com- 
pared to large spiral galaxies (e.g. Israel 2005). The possibility 
that molecular gas may be present in large quantities in such 



regions has been p ostulated ( see, e .g., IPfenniger et ail 1 1994 
CombeJUfflJ iKalherla et . alJbfXXl and references therein). 
Previous authors have concentrated on clumps with extents 
of the order of 10AU and masses very much smaller than a 
solar mass each. Such clumps would have a covering factor 
(or beam filling factor) that is much smaller than unity. Their 
existence might have a range of observational consequences. 
For instance, they may be responsible for "extreme scattering 
events", causing large variation s in ra dio signals from more dis- 
tant sources dWalker&Wardleiri 998). We focus on lower den- 
sity, larger objects having a covering factor near unity. 

Of course, no survey will lead to the detection of all neutral 
atomic or molecular gas in a galaxy. Sensitivity limits lead to 
mapped distributions having boundaries, even though the real 
gas distributions may be more extended. The purpose of the 
present investigation is the determination of parameters char- 
acterizing the gas and environment that could result in signifi- 
cant hidden quantities of gas, with column densities compara- 
ble to those often found by surveys to be associated with star 
formation. Such gas could provide a reservoir for eventual star 
formation which would in turn create conditions leading to the 
gas being more easily detected. In many cases, the observed 
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regions will be sufficiently distant that individual clouds are 
not observed. We have assumed, in effect, that all clouds have 
the same total column density and that, on average, any line of 
sight will pass through one cloud. 

Therefore, we consider the possible existence of molecular 
clouds that contain a significant column density of H2 whilst 
producing sufficiently little emission that they are undetectable 
or at least ignored by those conducting most surveys. We pro- 
duce theoretical models of the chemistry and emission arising 
in molecular clouds and conduct a parameter space search for 
physical conditions that might give rise to undetectable quan- 
tities of molecular gas. In Sect. |3] we describe the photon- 
dominated region (PDR) code used to calculate the chemical 
and emission properties of model clouds and outline the range 
of parameter space examined. In Sect. [3] we review the sen- 
sitivity limits of current telescopes appropriate to the detec- 
tion of emission (in the IR/sub-mm/mm/radio) from molecu- 
lar clouds and specify the constraints we impose for a cloud to 
be considered undetectable. The results of our parameter space 
search are presented in Sect.|4]and we summarize our findings 
in Sect.|5] 

2. Cloud model 

We u se the ucl_pdr time-dependent PDR code (see lBell et alJ 
120061 for a detailed description) to calculate the chemistry, 
thermal balance and emission strength of fine structure and 
molecular lines within a model cloud, for a given set of physi- 
cal parameters. 

A PDR is modelled as a one-dimensional semi-infinite slab, 
illuminated from one side by the standard Draine radiation field 
JPrainelll978l) scaled by a factor \- The chemistry and ther- 
mal balance are calculated self-consistently at each depth point 
into the slab and at each time-step, producing chemical abun- 
dances, emission line strengths and gas temperatures as func- 
tions of depth and time. Constant thermal pressure is assumed, 
and the density is calculated self-consistently with the tempera- 
ture at each depth- and time-step. The turbulent contribution to 
the pressure is neglected. Usually the temperature varies by at 
most a factor of a few in most models, and the assumption of a 
constant density would lead to similar results. Self-gravity may 
be important in real clouds, and lead to some density variation. 
However, magnetic forces may also be important and counter- 
act the gravity to a large extent. Given the uncertainities in the 
real density distributions, we consider that we are justified in 
making the isobaric assumption. 

We adopt a reduced chemical network containing 32 
species and over 300 reactions, including ion-molecule, pho- 
toionization and photodissociation reactions. Freeze-out of 
atoms and molecules onto grains is neglected. The gas-phase 
elemental abundances adopted are listed in Tableland are as- 
sumed to scale linearly with metallicity, Z (see below). The 
species included in the network are listed in Table |2 All frac- 
tional abundances in this paper are quoted relative to the hydro- 
gen nuclei number density [«h = "(H) + 2n(H?)"|. The reaction 
rates contained in the UMIST99 database jLe Teuff et aljE oOO) 
are adopted, with some modifica tions introduced a s part of a 
recent PDR benchmarking effort (Rollig et al. 2006). 



Table 1. Elemental abundances used in the reduced chemical 
network (relative to the hydrogen nuclei number density, «h)- 



He 


7.50xl(T 2 


C 


1.42x10-* 


O 


3.19xl(T 4 


Mg 


5.12xl(T 6 



Table 2. Species included in the reduced chemical network. 



H, H + , H 2 , Ht, H+, He, He + 

C, C + , CH, CH + , CH 2 , CH+, CH 3 , 

CH+, CH 4 , CH+, CHj, CO, CO + , HCO+ 

O, + , OH, OH + , H 2 0, H,0 + , H 3 + , 2 , O+ 

Mg, Mg + , e~ 



Variation in the metallicity (Z/Z ) is taken into account 
through the assumption that the dust and PAH number densi- 
ties, their photoelectric heating rates, the formation rate of H2 
on grains and the elemental metal abundances all scale linearly 
with Z. 

The ucl_pdr code calculates the integrated intensities of the 
[Oi] 63 ^im and 146 \im, [Ci] 610 ^im and 370 ^im, and [Cn] 
158 \xm fine structure lines and all CO rotational lines up to 
J — 11 — » 10. Theoretical antenna temperatures are calculated 
on the assumption that emission from the surface of the cloud is 
distributed equally over one hemisphere {2n steradians) and are 
integrated over the linewidth to give integrated line intensities 
(in Kkm s _1 ). 

2.1. Model parameter space 

Two sets of initial conditions are assumed for the chemical 
abundances in the models: purely atomic (all hydrogen in H° 
and all carbon in C°), or purely molecular (all hydrogen in H2 
and all carbon in CO). These represent the two extremes for 
the initial state of the cloud, forming either from tenuous ma- 
terial or from dense cores that have subsequently relaxed to a 
diffuse state, therefore bracketing the range of possible starting 
conditions. We envisage that the clouds that are initially molec- 
ular have originated from much denser material. Consequently, 
the outer layers of atomic hydrogen around them are narrow. 
Our PDR model computes the H/H2 transition accurately as 
these objects evolve. The abundances, temperature and emis- 
sion properties are calculated at each time-step for a total pe- 
riod of 1 Gyr in each cloud model. Models are calculated for 
the entire region of parameter space with both sets of initial 
abundances. 

Some correlation between the radiation field and the cosmic 
ray ionization rate would be expected for extragalactic obser- 
vations where a beam samples the average interstellar UV field 
and cosmic ray flux incident on a cloud, as the high-mass stars 
that emit the bulk of the far ultraviolet radiation also become 
supernovae that create remnants in which cosmic rays are ac- 
celerated. Therefore, the cosmic ray ionization rate per hydro- 
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Table 3. The range of values covered in the parameter space. 10 4 



Parameter Values taken for each parameter 



t 


10\ 10 2 , 10 8 , 10" 


yr 


z 


10°, 10 _1 , 10~ 2 , 10~ 3 


z G 


P/k 


10', 10 2 , 10 3 , 10 4 


cm" 3 K 


X 


io~ 3 , io~ 4 , io~ 5 


Draine 


Vturb 


0.05, 0.1,0.3, 0.5 


km s 



gen molecule (f) is assumed to scale linearly with the incident 
radiation field strength (^-Draine) as: 



^=1.3xlO- I7 ;r K 1 ] 



(D 



We choose a large region of parameter space but restrict 
attention to low radiation intensities (<1(T 3 Draine). A lower 
limit to the radiation field is given by the metagalactic ra- 
diati on field, which has been considered extensively recently 
by [S ternberg et alJ The spectral shapes of the standard 

Draine ( 1978) and the metagalactic radiation fields differ, but 
a comparison at lOOnm implies that the metagalactic field is 
the weaker by four orders of magnitude. Of course, that field 
varies from point to point depending on the distances to the 
nearest bright galaxies. We consider a range of four orders of 
magnitude in metallicity and gas pressure. Table[3]lists the in- 
dividual parameter values considered. Since regions where the 
radiation background is low compared to the solar neighbour- 
hood are far from stars which provide energy to stir their en- 
vironments, such regions are likely to be fairly quiescent. The 
level of turbulence in the clouds considered here is expected 
to be significantly lower than is typic ally found in standar d 
nearby Galactic giant molecular clouds (Willia ms et alll 995). 
We therefore include models with microturbulent velocities in 
the range 0.05-0.5 km s _I . The assumed value of this velocity 
is important for the heating and cooling of the clouds. It is also 
important in our considerations of detectability. We note that an 
ensemble of clouds possessing a cloud-cloud velocity disper- 
sion would produce emissi on lines that are broa der than those 
from the individual clouds; Wolf ire et al J (11993) included such 
a cloud-cloud velocity dispersion in their models. We have not, 
in part because regions where few supernovae occur will have 
small cloud-cloud velocity dispersions, and also in part to re- 
strict the range of parameter space to a manageable level. Our 
neglect of a cloud-cloud velocity dispersion leads to a more re- 
stricted region of parameter space corresponding to molecular 
dark matter that is difficult to detect. 

In total, we constructed 384 time-dependent cloud models 
to cover the specified parameter space. The next section con- 
tains the criteria used to determine what part of parameter space 
is associated with molecular dark matter. 

3. Emission constraints 



10 3 



^ 10 2 



10' 




10° 



10 4 



10 3 



^ io 2 - 
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Fig. 1. Contour plots of the maximum H2 column density 
produced by initially molecular cloud models with v tur b = 
0.05 km s that meet the criteria of undetectable emission (see 
Sect. Contour levels range from iNfa, < 10 21 cm -2 (white), 
through 10 21 , 10 22 and 10 23 ctrT 2 , to N Hl > 10 24 cirT 2 (black). 
The three plots show the pressure-metallicity plane of the pa- 
rameter space considered, for incident radiation field strengths 
of 10~ 3 , 10~ 4 and 10~ 5 Draine, respectively. 



18 cirT 2 



The distributions of neutral atomic and molecular clouds are 

studied predominantly with H 1 2 1 cm, [0 1] , [C 1] and [C 11] fine of neutral atomic hydrogen column densities of 1 . 6 x 1 

structure and CO rotational lines, as well as dust continuum for a velocity disp ersion of lOkms' 1 and a comparable veloc- 

radiation. 21cm surveys can lead to the six sigma detection itv resolution (e.g. lGiovaneHi et al .120051) . However, published 
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Table 4. Molecular models with v tur b = 0.05 km s~' that satisfy the criteria of undetectable emission and A^h, > 10 21 cirT 2 . 
The value of Nn 2 listed for each model is the peak H2 column density, reached at time ?high ■ The period of time during which 
A^h 2 ^ 10 21 cm -2 in each model is defined by f m ; n and f max . 



z 


y 
A 


P/k 




'high 


/ 

' mni 


/ 

' max 


(Zo) 


(Draine) 


(cm- 3 K) 


(cm- 2 ) 


(yr) 


(yr) 


(yr) 


10° 


lO" 3 


10 3 


2.9xl0 22 


7xl0 6 


0x10" 


lxl0 lJ 


10° 


10" 5 


10 4 


2.9xl0 22 


3xl0 7 


0x10° 


lxlO 9 


i -' 


1(T 3 


10 1 


2.4xl0 21 


lxlO 7 


6xl0 6 


lxlO 7 


10 -1 


10~ 3 


10 2 


2.2xl0 21 


4xl0 7 


7xl0 6 


lxlO 8 


irr 1 


10 4 


10' 


2.3xl0 21 


6x10 s 


8xl0 7 


8xl0 8 


10 _1 


10" 5 


10' 


6.0xl0 21 


3xl0 7 


0x10° 


lxlO 9 


10" 1 


10" 5 


10 2 


6.3xl0 21 


4xl0 3 


0x10° 


lxlO 9 


10" 1 


10" 5 


10 3 


1.2xl0 22 


lxlO 6 


0x10° 


lxlO 9 


10- 1 


10" 5 


10 4 


2.8xl0 23 


4xl0 7 


0x10° 


lxlO 9 


10 2 


10- 3 


10 1 


J.JA1U 




^xlf) 6 

JA1U 


lxlO 7 


10 2 


lO" 3 


10 2 


4.7X10 21 


8xl0 6 


4xl0 6 


lxlO 7 


lO" 2 


10" 3 


10 3 


1.6xl0 21 


5xl0 7 


2xl0 7 


lxlO 8 


10~ 2 


10- 4 


10 1 


5.5xl0 21 


9xl0 7 


0x10° 


lxlO 8 


10" 2 


10" 4 


10 2 


2.8xl0 21 


5xl0 8 


0x10° 


lxlO 9 




10- 4 


10 3 


l.lxlO 21 


lxlO 9 


0x10° 


lxlO 9 


lO" 2 


10 4 


10 4 


1.4xl0 21 


5xl0 3 


0x10° 


lxlO 7 


lO 2 


10 5 


10' 


2.6xl0 22 


0x10° 


0x10° 


lxlO 9 


lO" 2 


10"5 


10 2 


2.4xl0 22 


3xl0 5 


0x10° 


lxlO 9 


lO" 2 


10" 5 


10 3 


2.8xl0 22 


0x10° 


0x10° 


lxlO 9 


lO" 2 


10" 5 


10 4 


1.2xl0 24 


0x10° 


0x10° 


lxlO 9 


10-3 


10-3 


10 1 


1.2xl0 22 


9xl0 6 


2xl0 5 


lxlO 7 


in 3 


10 3 


10 2 


1 v 1 n 2 1 


lxlO 7 


2xl0 6 


lxlO 7 


1(T 3 


10- 3 


10 3 


7.9xl0 21 


9xl0 6 


3xl0 6 


lxlO 7 


10- 3 


10- 3 


10 4 


3.5xl0 21 


2xl0 7 


3xl0 6 


7xl0 7 


10- 3 


10- 4 


10' 


l.OxlO 22 


8xl0 7 


0x10° 


1x10 s 


10- 3 


10- 4 


10 2 


l.OxlO 22 


9xl0 7 


0x10° 


lxlO 8 


lO 3 


10- 4 


10 3 


9.1xl0 21 


2x10 s 


0x10° 


9xl0 8 


10- 3 


10- 4 


10 4 


3.3xl0 21 


1x10" 


0x10° 


lxlO 9 


lO 3 


10- 5 


10' 


2.1xl0 22 


9xl0 8 


0x10° 


lxlO 9 


10- 3 


10- 5 


10 2 


2.0xl0 22 


9xl0 4 


0x10° 


lxlO 9 


10- 3 


10- 5 


10 3 


2.0xl0 22 


2xl0 4 


0x10° 


lxlO 9 


10- 3 


10- 5 


10 4 


3.2xl0 22 


2xl0 5 


0x10° 


lxlO 9 



column density distributions of H 1 in other galaxies usua lly do 
not go to such low column densities. For instance. [Walter et alJ 
(2005J) take a column density of 10 20 ctrT 2 for the lowest con- 
tour in fi gures, and the lowes t column density shown in figures 
bv iNoordermeer et alJ J2005I) is roughly 10 19 cirT 2 . We shall 
use 10 19 crrT 2 as the minimum column density of neutral hy- 
drogen to be considered of observational significance. Though 
we have cited observations that were sensitive to Hi column 
densities roughly as low as 10 I9 cm -2 , we realize that this is a 
rather small value. We have chosen to define such a low limit- 
ing value for the H 1 column density because the small velocity 
dispersions that we have assumed would make H 1 more readily 
detectable. 

Ground-based observations of the lower CO rotational 
emission lines have a typical RMS noise limit of ~lmK; 
we therefore adopt a minimum level for the detectable in- 
tegrated line intensity of 0.01 Kkms~' (assuming a veloc- 
ity resolution of lOkms" 1 ), suitable for dee p searches (e.g. 
Iwilson & Combesll 19981 ISmoker et alJl200ol) . The noise lim- 



its on [Ci] 610 |xm ground-based observations are higher, but 
for convenience we also impose a constraint of 0.01 Kkms -1 
for the minimum detectable integrated line intensity. 

The [Oi] and [Cn] fine structure lines are only observ- 
able from space or the upper atmosphere. Taking typical RMS 
noise limits for spectra obtained with the Long Wavelength 
Spectrometer (LWS) on the Infrared Space Observatory (ISO), 
a minimum integrated flux of ~10~ 17 WrrT 2 is reasonable (as- 
suming a resolution of 0.29 [im for the [0 1] 63 \xm line and 
0.6 (J.m for the [C 11] 158 (J.m line). This corresponds to a surface 
brightness o f ~10~ 7 erg s' 1 cm' 2 sr -1 (assuming Ob ~ 10~ 7 sr 
for the LWS: lLlovdl2 003). We adopt this as the detection limit 
for the two dominant fine structure lines. 

Since the dust in these clouds is cold (typically <5 K) and 
depleted at the lower metallicities, its emission will be faint and 
therefore unlikely to be detected above the background level. 

To summarize, the limits for non-detection that we im- 
pose on the cloud models are: A^hi < 10 19 crrT 2 , Wco < 
0.01 Kkms~' in all rotational lines, W < 0.01 Kkms 1 for 
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10 ' HP 10" 5 10 ; 10 ; 10" 5 

X (Draine) % (Draine) 

Fig. 2. Contour plots of the maximum H2 column density produced by initially molecular cloud models with Vturb = 0.05 km s 
that meet the criteria of undetectable emission (see Sect. [5}. Contour levels range from Nn 2 < 10 21 crrT 2 (white), through 10 21 , 
10 22 and A^h, > 10 23 cm -2 (dark-grey). The four plots show the pressure-radiation plane of the parameter space considered, for 
metallicities of 1, 10 _1 , 10~ 2 and 10~ 3 times solar (from left to right, top to bottom). 



the [Ci] 610 |xm line, / < 10~ 7 erg s _1 cirr 2 sr~' for the [Oi] 
63 [im and [C 11] 158 \xm lines. Any cloud model with emission 
below all these limits simultaneously is considered to be effec- 
tively undetectable with current technology. We now describe 
the results of our parameter space search, presenting models 
that meet these criteria. 

4. Model results 

The column density through the cloud is restricted in each 
model such that the total emission in each line is below the cor- 
responding limit given above. For each set of physical param- 
eters, the resulting cloud model therefore represents the maxi- 
mum amount of molecular gas that can be hidden under those 
conditions. We then select cloud models that contain signifi- 
cant quantities of H2, specifically > 10 21 crrr 2 , at some 
point during their 1 Gyr evolution. Such a column density is of 
the same order of magnitude as the hydrogen nuclei column 
densities associ ated with parts of l arge spiral galaxies in which 
stars form (e.g.lNoordermee ret alJl2005l) . This subset of mod- 
els defines the region of our parameter space where significant 
amounts of molecular gas can reside without being detected in 



emission. Tables |4j{7]list the subset of cloud models with ini- 
tially molecular gas that meet our criteria for non-detectable 
emission whilst containing a significant amount of H2. Table|S] 
lists the subset of initially atomic cloud models that meet the 
same conditions. The tables also show the period of time for 
which this is the case in each model, starting at f m i„ and ending 
at ? m ax, as well as the time, fygh, at which the peak value of A^h, 
is reached. In total, 1 1 1 of the 384 model clouds obtain H2 col- 
umn densities >10 21 crrr 2 for some period of their evolution. 
The trends describing the change in total H2 column density 
with parameter variation are discussed in the next two sections. 

4.1. Initially molecular clouds 

For models for which the gas is initially molecular, the peak 
H2 column density (>10 24 cirT 2 ) occurs in the cloud with the 
lowest turbulent velocity, subject to the lowest incident radia- 
tion field strength, with low metallicity and with the highest 
pressure (0.05 km s~\ 10 -5 Draine, 1O~ 2 Z , 10 4 cirT 3 K; see 
Tableland Fig.^. For the lowest radiation field strength, in 
the models with low turbulent velocities (0.05 or 0.1 kms~'), 
the H2 column density increases as the pressure rises and it 
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Table 5. Molecular models with Vturb = O.lkms -1 that satisfy the criteria of undetectable emission and Nn 2 > 10 21 cm -2 . 
The value of Nn 2 listed for each model is the peak H2 column density, reached at time ?high ■ The period of time during which 
Nft 2 > 10 21 cirT 2 in each model is defined by f mm and f max . 



z 


X 


P/k 


N S2 




tmm 






(Draine) 


(cm -3 K) 


(cm -2 ) 


(vr) 


(vr) 


(vr) 


10° 


10- 5 


10 4 


2.9xl0 22 


2xl0 6 


OxlO u 


lxlO lJ 


10" 1 


10- 3 


10 1 


1.4xl0 21 


9xl0 6 


8xl0 6 


IxlO 7 


10- 1 


10" 3 


10 2 


2.6xl0 21 


4xl0 7 


IxlO 7 


IxlO 8 


10" 1 


10- 4 


10' 


2.5xl0 21 


5xl0 8 


IxlO 8 


IxlO 9 


10 -1 


10" 5 


10' 


2.7xl0 21 


6xl0 5 


0x10° 


IxlO 9 


10" 1 


10 5 


10 2 


2.7xl0 21 


0x10° 


0x10° 


IxlO 9 


10- 1 


10" 5 


10 3 


3.1xl0 21 


2xl0 5 


0x10° 


IxlO 9 


10" 1 


10" 5 


10 4 


7.0xl0 22 


5xl0 5 


0x10° 


IxlO 9 


io~ 2 


10- 3 


10 1 


5.0xl0 21 


9xl0 6 


4xl0 6 


IxlO 7 


10" 2 


10" 3 


10 2 


4.0xl0 21 


8xl0 6 


4xl0 6 


IxlO 7 


10" 2 


10~ 3 


10 3 


1.7xl0 21 


7xl0 7 


2xl0 7 


9xl0 7 


10" 2 


10" 4 


10 1 


4.5xl0 21 


8xl0 7 


4xl0 7 


IxlO 8 


10" 2 


10" 4 


10 2 


3.0xl0 21 


4xl0 8 


6xl0 7 


IxlO 9 


10" 2 


10" 4 


10 3 


l.OxlO 21 


9xl0 8 


8xl0 8 


9xl0 8 


10-2 


10" 5 


10' 


3.0xl0 21 


9xl0 3 


0x10° 


IxlO 9 


10" 2 


1(T 5 


10 2 


2.9xl0 21 


0x10° 


0x10° 


IxlO 9 


10~ 2 


10 5 


10 3 


3 IxlO 21 


2xl0 6 


0x10° 


IxlO 9 


10" 2 




10 4 


4.6xl0 21 


7xl0 5 


0x10° 


IxlO 9 




1U 


10 1 


1.2xl0 22 


IxlO 7 




IxlO 7 


10 3 


10- 3 

1U 


10 2 


l.lxlO 22 


IxlO 7 


2xl0 6 


IxlO 7 


10- 3 


10- 3 


10 3 


6.8xl0 21 


9xl0 6 


3xl0 6 


IxlO 7 


10- 3 


10- 3 


10 4 


4.0xl0 21 


3xl0 7 


5xl0 6 


8xl0 7 


10- 3 


10- 4 


10' 


l.lxlO 22 


1x10 s 


2xl0 6 


IxlO 8 


10- 3 


10- 4 


10 2 


l.OxlO 22 


9xl0 7 


8xl0 6 


IxlO 8 


10- 3 


10- 4 


10 3 


7.8X10 21 


2x10 s 


2xl0 7 


IxlO 9 


10- 3 


10- 4 


10 4 


3.0xl0 21 


9xl0 8 


5xl0 7 


IxlO 9 


10- 3 


10 s 


10' 


l.lxlO 22 


9xl0 8 


0x10° 


IxlO 9 


10- 3 


10- 5 


10 2 


8.1xl0 21 


IxlO 9 


0x10° 


IxlO 9 


10- 3 


10- 5 


10 3 


3.1xl0 21 


1x10" 


0x10° 


IxlO 9 


10- 3 


10- 5 


10 4 


1.7xl0 21 


4xl0 5 


0x10° 


IxlO 9 



is the intermediate metallicities (10 _1 -10 _2 Z o ) that produce 
the highest column densities of dark material. For models for 
which the radiation field does not have its lowest value, the H2 
column density increases as the pressure and metallicity drop 
(see Fig.^. 

At the lowest metallicity considered (10~ 3 Z o ), the H2 col- 
umn density is >10 21 cirT 2 for almost all combinations of tur- 
bulent velocity, radiation field strength and pressure (see the 
bottom right panel of Fig.|5]for the v tul b = 0.05 km s~' case). 
The only exceptions are the three cloud models with the lowest 
incident radiation and highest pressure at turbulent velocities of 
0.3 and 0.5 km s~' (see Tables[6]and[7}. Overall, increasing tur- 
bulent velocity within the cloud serves to reduce the H2 column 
density that can be hidden, although some of the low pressure 
models show a slight increase in A^h, at higher turbulent veloc- 
ities. 

The period of time for which Nh 2 ^ 10 21 cm 2 in each 
model (i.e. from to t mm ) is governed by the restrictions 
placed on the CO and Hi emission. For cloud models with 
fmin > 0, the 2.6 mm CO line emission is initially above the 
limit given in Sect. |3] and therefore considered detectable, but 



falls below the limit once f m ; n is reached. In models with 
fmax < 1 Gyr, the maximum time considered in the calcula- 
tions, the 21 cm Hi emission rises above its detection limit at 
late times, reducing the period for which the cloud remains un- 
detectable. This is generally the case for models with higher 
radiation field strengths. 



An analysis of the thermal balance in the models shows that 
the dominant cooling processes are emission in the lower rota- 
tional transitions of CO and the fine structure lines of neutral 
carbon, with some contribution from collisions with the cooler 
dust grains in the high pressure models. The heating of the gas 
at the cloud surface is dominated by H2 formation and pho- 
todissociation, whilst turbulence and cosmic ray ionization be- 
come important heating mechanisms deeper into the cloud. In 
models with higher radiation field strengths, PAH photoelec- 
tric heating plays a role near the cloud surface. Heating due to 
turbulent dissipation becomes the dominant heating process at 
higher turbulent velocities. 
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Table 6. Molecular models with v tur b = 0.3kms _1 that satisfy the criteria of undetectable emission and A 7 ^ > 10 21 cm -2 . 
The value of A 7 ^ listed for each model is the peak H2 column density, reached at time ?high ■ The period of time during which 
Nn 2 > 10 21 cirT 2 in each model is defined by f min and f max . 



z 




P/k 




b nign 


/ ■ 

' nun 


/ 

' max 


(Zn.) 


(Draine) 


(cm -3 K) 


fcrrr 2 ) 


(vrl 




(vt) 


io-' 


10- 3 


10 1 


2.0xl0 21 


5xl0 6 


5xl0 6 


lxlO 7 


10- 1 


10- 3 


10 2 


2.3xl0 21 


1x10 s 


7xl0 6 


lxlO 8 


10" 1 


10 4 


10' 


3.4xl0 21 


9xl0 8 


4xl0 7 


lxlO 9 


10 2 


10- 3 


10 1 


5.6xl0 21 


9xl0 6 


3xl0 6 


lxlO 7 


10" 2 


10" 3 


10 2 


4.4xl0 21 


9xl0 6 


3xl0 6 


lxlO 7 


10 -2 


10" 3 


10 3 


1.5xl0 21 


3xl0 7 


5xl0 6 


lxlO 8 


10~ 2 


10" 4 


10 1 


5.1xl0 21 


7xl0 7 


3xl0 7 


lxlO 8 


10" 2 


10" 4 


10 2 


5.5xl0 21 


2xl0 8 


4xl0 7 


lxlO 9 


IO" 2 


10" 5 


10 1 


5.3xl0 21 


lxlO 9 


4xl0 8 


lxlO 9 


10" 2 


10 5 


10 2 


1.7xl0 21 


lxlO 9 


5xl0 8 


lxlO 9 


10- 3 


10- 3 


IO 1 


1.2xl0 22 




0x10° 


lxlO 7 






1 u 


1 .ZX 1U 


1 N/ 1 CP 
1 X 1U 


On/ 1 f»5 
ZX 1U 


1 N/ 1 cP 
1 X 1U 


10- 3 


10- 3 


10 3 


6.0xl0 21 


9xl0 6 


3xl0 6 


lxlO 7 


10- 3 


10- 3 


10 4 


1.9xl0 21 


7xl0 6 


5xl0 6 


lxlO 8 


10- 3 


10- 4 


10' 


1.3xl0 22 


lxlO 8 


0x10° 


lxlO 8 


10- 3 


10- 4 


10 2 


9.3xl0 21 


9xl0 7 


2xl0 6 


1x10 s 


10- 3 


10- 4 


10 3 


3.9xl0 21 


7xl0 7 


3xl0 7 


lxlO 9 


10- 3 


io~ 4 


10 4 


2.2X10 21 


6x10 s 


8xl0 7 


lxlO 9 


10- 3 


io 5 


10' 


1.2xl0 22 


1x10" 


0x10° 


lxlO 9 


10- 3 


10- 5 


10 2 


l.OxlO 22 


1x10" 


2xl0 7 


lxlO 9 


10- 3 


10 5 


10 3 


3.5xl0 21 


6xl0 8 


3xl0 8 


lxlO 9 



4.2. Initially atomic clouds 



From Table [S] it can be seen that there are significantly fewer 
models for which the gas is initially atomic that can produce 
appreciable quantities of hidden molecular gas. This is mainly 
due to the relatively long time-scales required for molecular 
hydrogen formation under the physical conditions being con- 
sidered here. However, despite this limitation, there still exists 
a well-defined region of parameter space within which atomic 
gas can evolve to form large column densities of molecular gas 
whilst remaining undetectable in emission by the usual tracers. 

In all the models listed in Table [S] the radiation field 
strength is the lowest considered in our parameter space 
(l(T 5 Draine) and the pressures are among the highest (10 3 
or 10 4 crrT 3 K). There is some spread in metallicity, but there 
are no clouds with turbulent velocities above O.lkms -1 that 
have >10 21 crrT 2 of H2. The peak H2 column density is 
7.4xl0 22 cm~ 2 , some 20 times smaller than the largest column 
density obtained by clouds evolving from molecular gas. None 
of the initially atomic clouds reach A 7 ^ > 10 21 cirT 2 until they 
have existed for 20 Myr or longer; several of the cloud mod- 
els do not obtain significant H2 column densities until 1 Gyr 
has passed. The heating and cooling of the gas is governed by 
the processes discussed in the previous section, although the 
heating due to the formation of H2 is more significant in these 
models, since the hydrogen is initially atomic. 



5. Conclusion 

We have conducted a large study of the parameter space asso- 
ciated with regions where the radiation field is weak relative to 
that in the solar vicinity to search for conditions that give rise to 
significant quantities of molecular gas with emission strengths 
that are comparable to the lower limits considered in standard 
surveys. Having explored parameter space to identify the time- 
dependent cloud models that satisfy these limits, we selected 
those models in which a significant column density of molecu- 
lar hydrogen, A 7 ^ > 10 21 cm' 2 was reached. 

We find that there is a large region within the parame- 
ter space that meets these requirements, resulting in clouds 
that might contain a significant mass of molecular gas whilst 
remaining effectively undetectable or at least not particu- 
larly noticeable in surveys. This region is characterized by 
low turbulent velocities (<0.1kms _1 ), low radiation fields 
(<10~ 4 Draine), normal to low metallicities (~1-1O~ 2 Z ) and 
intermediate pressures (~10 3 -10 4 cirT 3 K). Clouds forming 
from molecular gas under these conditions are capable of main- 
taining large amounts of H2 without becoming detectable in 
emission. Initially atomic gas can also form significant column 
densities of molecular material with undetectable emission, but 
the time-scales required are much longer (>10 7 yr). 

Sternberg et al. (2002) conducted a study of PDRs in which 
they assumed the present metagalactic radiation field to be in- 
cident on clouds in work on minihaloes composed of cold dark 
matter of the variety considered in many cosmological studies 
and on compact high velocity H 1 clouds. As mentioned above, 
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Table 7. Molecular models with Vturb = 0.5kms _I that satisfy the criteria of undetectable emission and Nn 2 > 10 21 cm -2 . 
The value of Njj 2 listed for each model is the peak H2 column density, reached at time ?high ■ The period of time during which 
A^h 2 > 10 21 cirT 2 in each model is defined by f min and f max . 



z 




P/k 




'nign 


/ 

' nun 


/ 

' max 


(Zn.) 


(Draine) 


(cm -3 K) 


(crrr 2 ) 


(vr) 




(vt) 


10 1 


10- 3 


10' 


1.5xl0 21 


lxlO 7 


9xl0 6 


lxlO 7 


10- 1 


10- 3 


10 2 


2.8xl0 21 


6xl0 7 


8xl0 6 


lxlO 8 


lO" 1 


10~ 4 


10' 


3.9xl0 21 


6xl0 8 


9xl0 7 


lxlO 9 


10 2 


10- 3 


10 1 


4.0xl0 21 


lxlO 7 


4xl0 6 


lxlO 7 


10" 2 


10" 3 


10 2 


5.3xl0 21 


lxlO 7 


4xl0 6 


lxlO 7 


10~ 2 


10" 3 


10 3 


1.3xl0 21 


2xl0 7 


6xl0 6 


lxlO 8 


10~ 2 


10" 4 


10 1 


3.8xl0 21 


lxlO 8 


4xl0 7 


lxlO 8 


10" 2 


10" 4 


10 2 


5.7xl0 21 


4xl0 8 


5xl0 7 


lxlO 9 


lO" 2 


10" 5 


10 1 


3.4xl0 21 


lxlO 9 


4xl0 8 


lxlO 9 


10" 2 


10 5 


10 2 


2.2xl0 21 


lxlO 9 


6xl0 8 


lxlO 9 


io- J 


10- 3 


10 1 


9.4xl0 21 


lxlO 7 


0x10" 


lxlO 7 


10- 3 


10- 3 


10 2 


8.0xl0 21 


lxlO 7 


2xl0 5 


lxlO 7 


10- 3 


10- 3 


10 3 


5.1xl0 21 


9xl0 6 


3xl0 6 


lxlO 7 


10- 3 


10- 3 


10 4 


l.lxlO 21 


5xl0 6 


5xl0 6 


9xl0 7 


10- 3 


10- 4 


10' 


l.lxlO 22 


lxlO 8 


0x10° 


lxlO 8 


10- 3 


10- 4 


10 2 


8.4xl0 21 


1x10 s 


2xl0 6 


1x10 s 


10- 3 


10- 4 


10 3 


3.9xl0 21 


8xl0 7 


2xl0 7 


lxlO 9 


10- 3 


10- 5 


10' 


l.OxlO 22 


lxlO 9 


0x10° 


lxlO 9 


10- 3 


10- 5 


10 2 


8.3xl0 21 


1x10" 


8xl0 6 


lxlO 9 


10- 3 


10- 5 


10 3 


3.3xl0 21 


9xl0 8 


3xl0 8 


lxlO 9 



Table 8. All atomic models that satisfy the criteria of undetectable emission and A^h 2 > 10 21 cm -2 . The value of A^h, listed for 
each model is the peak H2 column density, reached at time fhigh- The period of time during which A^h 2 > 10 21 cm -2 in each model 
is defined by t m [ n and f max . 



Vturb 


Z 


X 


P/k 


Nn 2 


fhish 




►max 


(kms" 1 ) 


(Z Q ) 


(Draine) 


(cm- 3 K) 


(cm~ 2 ) 


(yr) 


(yr) 


(yr) 


0.05 


10° 


10- 5 


10 3 


1.8xl0 21 


2xl0 8 


2xl0 8 


1x10" 


0.05 


10° 


10- 5 


10 4 


2.9xl0 22 


2xl0 7 


2xl0 7 


1x10" 


0.05 


10- 1 


10- 5 


10 3 


2.1xl0 21 


1x10" 


9x10 s 


lxlO 9 


0.05 


10- 1 


10- 5 


10 4 


7.4xl0 22 


3xl0 8 


2xl0 8 


lxlO 9 


0.05 


10- 2 


10- 5 


10 4 


l.lxlO 21 


lxlO 9 


1x10" 


lxlO 9 


0.1 


10° 


10- 5 


10 4 


2.9xl0 22 


2xl0 7 


2xl0 7 


lxlO 9 


0.1 


10- 1 


lO" 5 


10 3 


1.5xl0 21 


lxlO 9 


1x10" 


lxlO 9 



the radiation field that they assumed is of comparable strength 
at lOOnm to fields for which we found large amounts of dark 
molecular gas to exist. Hence, there is some overlap between 
their work and ours. However, their emphasis is on Hi clouds 
detectable in 21 cm emission, whereas ours is on the possible 
existence of dark molecular gas with substantial enough col- 
umn densities to serve as reservoirs for potential star formation. 
They address observations of H 1 clouds having total hydrogen 
nuclei column densities of an order of magnitude and more 
below those with which we are concerned. Our interest is in 
whether regions of much higher column density may harbour 
sufficient molecular gas to become sites of stellar birth while 
they remain hidden in current surveys. We also wonder whether 
molecular dark matter existing in regions far away from bright 
stars may lead to future star formation in some galaxies oc- 



curring significantly beyond the currently detected bounds of 
those galaxies. 

Finally, we note that an H2 column density of 10 21 cm -2 
distributed evenly over a disc with a radius of 40kpc would 
correspond to a mass of 10 11 M Q , sufficient to be of dynamical 
importance in some galaxies. 

A means of placing limits on the product of the covering 
factor of the type of molecular dark matter we are consider- 
ing and the mean surface area per galaxy over which such 
dark matter is distributed would be through millimetre-wave 
absorptio n observations of other galaxies similar to those made 
bv iLiszt & Lucasl 1 1 998 ) to study diffuse clouds in our galaxy. 
Observations towards many background sources would almost 
certainly be required to set meaningful limits on that product, 
as a detection rate of about 4 optical QSO absorption systems 
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per unit redshift at a redshift of ~1 implies that intervening 
galaxies must have r adii of almost 100 kp c if they are randomly 
oriente d discs (e.g. ISareent et alJll979l) . IZwaan & Proch aska 
(2006) have argued that a redshift interval of Az ~ 100 has 
been covered in the CORALS survey of radio-selected quasars 
(Ellis on etai]l200ll) . whereas a Az range on the order of 3300 
near z = should be required to discover a high H2 column 
density absorber at the present epoch, if H2 has the sort of dis- 
tributions in galaxies that are i nferred from the BIMA SONG 
maps of CO in nearby galaxies ( Heife r et al.l2003l) . 
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